Scaling laws for the photo-ionisation cross section of two-electron atoms 
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The cross sections for single-electron photo-ionisation in two-electron atoms show fluctuations 
which decrease in amplitude when approaching the double-ionisation threshold. Based on semiclas- 
sical closed orbit theory, we show that the algebraic decay of the fluctuations can be characterised in 
terms of a threshold law a oc as E —* 0- with exponent fi obtained as a combination of stability 
exponents of the triple-collision singularity. It differs from Wannier's exponent dominating double 
ionisation processes. The details of the fluctuations are linked to a set of infinitely unstable classical 
orbits starting and ending in the non-regularisable triple collision. The findings are compared with 
quantum calculations for a model system, namely collinear helium. 
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PACS numbers: 32.80.Fb,03.65.Sq,05.45.Mt,05.45.-a 

Recent experimental progress has significantly im- 
proved the energy resolutions of highly excited two elec- 
tron states below [l|,|2| an d above [3[ the double ionisation 
threshold E — 0. Near the threshold, electron-electron 
correlation effects become dominant which are directly 
observable in total and partial cross sections, see 0, d, @| 
for recent reviews; an example is Wannier's celebrated 
threshold law for double ionisation Q confirmed experi- 
mentally in Q. For E < 0, two-electron atoms exhibit a 
rich resonance spectrum while the classical dynamics of 
this three-body Coulomb problem becomes chaotic. Ap- 
proaching the double ionisation threshold from below has 
thus proved challenging [f| and recent experimental and 
theoretical efforts still only reach single-ionisation thresh- 
olds In with ./V < 15 0,10. 

Semiclassical methods need 
to address the chaotic nature of the classical dynamics 
which is dominated by the complex folding patterns of 
the stable/unstable manifolds of the triple collision pjj |. 
Due to the high dimensionality of the system semiclas- 
sical applications have been restricted to subsets of the 
full spectrum and again small N values [f| . 

In this letter, we show that the fluctuations in the total 
cross section for single electron photo-ionisation below 
the three-particle breakup energy decays algebraically 
with an exponent determined by the triple collision sin- 
gularity different from Wannier's exponent. Writing the 
cross section in dipole approximation (ill . [l2| , we obtain 



a{E) = -4tt atujJ^{D<f>i\G{E)\D^i 



(1) 



where <pi is the wave function of the initial state and 
D = 7r • (ri + Y2) is the dipole operator with tt, the 
polarisation of the incoming photon and Tj , the position 
of electron j. Furthermore, G(E) is the Green function of 
the system at energy E = Ei + hui and a = e 2 /He. Note 
that we work in the infinite nucleus mass approximation, 
that is, the position of the nucleus is fixed at the origin. 

By expressing the Green function semiclassically in 
terms of classical trajectories [ll|, fluctuations in the 
cross section of hydrogen-like atoms in external fields 



have been analysed successfully using closed orbit theory 
(COT) [11IIT2I]. In the semiclassical limit, the support of 
the wave function fa shrinks to zero relative to the size 
of the system reducing the integration in JT]) to an evalu- 
ation of the Green function at the origin. This is strictly 
valid only for potentials sufficiently smooth at the ori- 
gin, corrections due to diffractive scattering at the central 
singularity give additional contributions often treated in 
quantum defect approximation [ID, [l4[ . The situation is 
different for two-electron atoms where the dynamics near 
the origin is dominated by the non-regularisable triple 
collision singularity. Closed orbit theory has been used 
to analyse experimental photoabsorbtion spectra of he- 
lium with and without external fields [l3l Il5t | for highly 
asymmetric states; accompanying theoretical considera- 
tions treat the system in a single electron approximation 
thus not considering the triple collision dynamics which 
becomes important for doubly-excited resonances. 

In the following, we will discuss a COT treatment of 
two-electron atoms explicitly including the triple colli- 
sion dynamics when approaching the limit E — > 0_. In- 
troducing the hyperradius R = (r 2 + r 2 ,) 1 / 2 , we fix Rq 
such that a surface E defined as R — Rq encloses the 
support of initial state fa. The surface E naturally leads 
to a partition of the configuration space into physically 
distinct regions. In particular, quantum contributions to 
([T]) from the inner region are insensitive to the total en- 
ergy. Contributions from the outer region test the full 
scale of the classically allowed region of size \E\~~ 1 and 
will be responsible for the resonance structures near the 
double ionisation threshold E = 0. Following Granger 
and Greene [l6| , we write the photo-ionisation cross sec- 
tion ((T|) in terms of local scattering matrices, that is, 
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Here, S 1 ^ is a core-region scattering matrix which maps 
amplitudes of waves coming in at E onto amplitudes of 
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the wave components going out at E; it thus contains 
all the information about the correlated two-electron dy- 
namics near the nucleus. Likewise, describes the wave 
dynamics of the two-electron wave function emanating 
from and returning to E and thus picks up long-range cor- 
relation in the exterior of E. Furthermore, d is the atomic 
dipole vector, d n {E) =< ^(E)\D<f>i > with ^^(E) be- 
ing the nth linearly-independent energy-normalised so- 
lution of the Schrodinger equation inside E with incom- 
ing wave boundary conditions at E [16]. This type of 
scattering formulation was independently developed in 
[ijj for general surfaces E, for a semiclassical formula- 
tion see 1181 . The series expansion ([3]) was exploited in 
(l3l LLJ, LL9| in order to include core-region scattering or 
quantum defect effects in COT. 

For the wave dynamics inside E, the double ionisation 
threshold E = is irrelevant and the core-region scatter- 
ing matrix as well as the dipole vector d vary smoothly 
across E — 0; they can be regarded as constant for ener- 
gies sufficiently close to the threshold. The information 
about the increasing number of overlapping resonances 
near the threshold is thus largely contained in . 

Semiclassical approximations for the quantities intro- 
duced above are valid in the outer-region R > Rq. The 
long-range scattering matrix, Sft , can thus be treated 
semiclassically while d and demand a full quantum 
treatment. The semiclassical representation of S' in po- 
sition space reads [ill E3, El] 

SHx,x',E) « (27rtfi)-' £ i 1 2|M ia |7 1 /2 e i*-*^ , (4) 
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where the sum is taken over all classical paths j con- 
necting points x and x 1 on E without crossing E; Sj(E) 
denotes the action of that path, Vj is the Maslov in- 
dex and / = 4 is the dimension of the system for 

— 1/2 

fixed angular momentum. Furthermore, |Mi2|j = 
\d 2 Sj(x, x^/dxdx'] 1 / 2 , where M12 is a (3 x 3) sub-matrix 
of the 6-dim. Monodromy matrix describing the lin- 
earised flow near a trajectory. Note that due to the strong 
instability of the classical dynamics near the triple col- 
lision, these matrix elements become singular for triple 
collision orbits (TCO) starting from or falling into the 
triple collision R = 0. It is thus important here that tra- 
jectories contributing to Q start at a fixed hyper-radius 
Rq > away from the triple collision (20j . 

Making use of the scaling properties of the classical 
dynamics, we introduce the transformation plj 

r = r/|£|; p = vl^P; S = S/y/\E\, L = l/y/\E\ 

where f , p corresponds to coordinates and momenta at 
fixed energy E = — 1 and L is the total angular momen- 
tum. Expressing L in scaled coordinates, we have L — > 
as E — > and can thus restrict the analysis to the three 
degrees of freedom subspace L — (for fixed L) [13]. 

In scaled coordinates, the inner-region shrinks accord- 
ing to Rq = |-E|i?o —> for E — > 0_, and the part of the 



dynamics contributing to in ((4]) is formed by trajec- 
tories starting and ending closer and closer to the triple 
collision R = as \E\ oc Rq — » 0. TCOs only occur in 
the so-called eZe space [l0( , a collinear subspace of the 
full three body dynamics where the two electrons are on 
opposite sides of the nucleus [21(. As Rq — > 0, only orbits 
coming close to the eZe space can start and return to E 
and they will do so in the vicinity of a closed triple col- 
lision orbit (CTCO) starting and ending exactly in the 
triple collision. The dynamics in the eZe space is rela- 
tively simple as it is conjectured to be fully chaotic with 
a complete binary symbolic dynamics. In particular, for 
every finite binary symbols string there is a CTCO, the 
shortest being the so-called Wannier orbit (WO) of sym- 
metric collinear electron dynamics. Furthermore TCOs 
escape from or approach the triple collision at R = 
always symmetrically along the T\ — ri axis in the eZe 
space, that is, along the WO [12] • This universality will 
be exploited below when treating the energy-dependence 
of Mia in 

Returning to the cross-section ©, we write a = o-q + 
Of 1, where we identify the smooth contribution o~q with 
the leading term in the series expansion ([3]). The main 
contribution to the fluctuating part of the cross section 
a fi is contained in S' which in semiclassical approxima- 
tion ([4]) can be expressed in terms of orbits returning to 

E once; multiple traversals of E represented by (S^S^) 
with k > 2 will give sub-leading contributions in the semi- 
classical limit E — > 0_ due to the unstable dynamics near 
the triple collision. Furthermore, swarms of trajectories 
starting on E and returning to E will do so close to the 
eZe subspace and thus in the neighbourhood of a CTCO 
with actions and amplitudes approaching those of the 
CTCO trajectory as Rq — > 0. The fluctuating part can 
thus in leading order be written in the form 

a fl (E)^n J2 ME)e uSi (5) 

CTCOj 

with 

ME) oc \M 12 (E)\J 1/2 = \E\ 9 ^\M 12 (R )\- 1/2 (6) 

and z = l/hy/\E\. In ([5|), the sum is taken over all 
CTCO's j starting and ending at E. Note that the sta- 
bility Mi2(i?o) in scaled coordinates depends on energy 
implicitly through the scaled radius Ro(E) = \E\Rq. 
As E — > 0_, M12 picks up additional contributions of 
parts of the CTCO closer and closer to the triple col- 
lision. Asymptotically, all CTCOs approach the triple 
collision along the WO and the contributions to M12 be- 
come orbit-independent. The i?-dependence for the Mon- 
odromy matrix of the WO can for small R be obtained 
analytically (23| leading to 

|M 12 (i?o)| oc |i? r 2A1+9/2 for R -> (7) 
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with exponent 



100Z-9 /4Z-9 



4Z- 1 



- 1 



(8) 



Thus, in unsealed coordinates, Mx% in (JSJ) diverges which 
is a direct consequence of the non-regularisability of the 
triple collision acting as an infinitely unstable point in 
phase space; details will be presented in [23[. The expo- 
nent fjj in ([5J) is universal for all CTCOs and consists of 
two components: /i e z e is related to the linearised dynam- 
ics in the eZe space and fi wr picks up contributions from 
two equivalent expanding degrees of freedom orthogonal 
to the eZe space in the so-called Wannier ridge (WR). 
The latter is the invariant subspace of symmetric elec- 
tron dynamics with |ri| = |7~2 1 at all times 21]. The fluc- 
tuations in the photo-ionisation cross section thus vanish 
in amplitude as E — > 0_ according to 



E 

CTCO, 



„iz.Sj 



(9) 



~^Aj depend only 



where the rescaled amplitudes a,j = \E\ 
weakly on E. These amplitudes contain contributions 
from the linearised dynamics along the orbit far from E 
as well as information about the inner quantum region 
R < Ro via the dipole vector d and the core-region scat- 
tering matrix SK Furthermore, multiple traversals of £ 
contained in (S^S^) with k > 2 in ([3]) approach the 
triple collision k times from the semiclassical side and 
will thus contribute at lower order with weights scaling 
at least as A^j ~ \E\ k ^ . The exponent /i in ([5]) is different 
from Wannicr's exponent fj, w with 



1 /100Z-9 



Hw - 4 



AZ- 1 



(10) 



One obtains, for example, fi = 1.30589... compared to 
fi w = 1.05589... for helium; the WR contributes to the 
decay for Z > 9/4 when fi wr is real. 

The exponent (i can be interpreted in terms of stabil- 
ity exponents of the triple collision singularity. Using an 
appropriate scaling of space and time by, for example, 
employing McGehee's technique 22], the dynamics near 
the singularity is dominated by two unstable fixed points 
in scaled phase space, the double escape point (DEP) and 
the triple collision point (TCP). In unsealed coordinates, 
these fixed points correspond to the WO at energy E = 0, 
that is, the DEP is the trajectory of symmetric double es- 
cape while the TCP corresponds to the symmetric triple 
collision and is the time reversed of the DEP. The triple 
collision itself can be mapped onto the classical dynam- 
ics at E — 0; likewise, i?o = l-^l-Ro acts as a parameter 
measuring the closeness to the E = manifold which 
contains the fixed points, see [13, H2] ■ Most classical tra- 
jectories emerging from E in the vicinity of the triple 
collision R = will lead to immediate ionisation of one 
electron carrying away a larger amount of kinetic energy. 



Only a fraction of orbits starting on £ near the WO will 
enter a chaotic scattering region; the WO at E < thus 
acts as an unstable direction of the DEP, U]j°, with a 
stability exponent . From there, they can return to 
the triple collision region and thus approach the surface 
£ again along the WO, that is, along the stable direction 
S^°. The transition from the DEP into the chaotic scat- 
tering region and from this scattering region to the TCP 
is limited by the least stable eigendirections of the fixed 
points in each of the invariant subspaces (eZe or WR) 
perpendicular to the WO. 

Trajectories leaving £ along the WO in the eZe space 
diverge from the WO along an unstable direction Uff e 
with a stability exponent A^? e . Competition of the insta- 
bility in Uf) with that in J7|) Ze determines the fraction 
A® Ze of orbits entering the chaotic scattering region. By 
using methods as in [l(| [24[ , one finds that 



-ttg 



(11) 



with exponent equal to Wannier's exponent (p~0|) . For the 
cross section (jTJ , information about the phase space re- 
gion returning from the chaotic scattering region to the 
surface £ along the WO is also needed. While approach- 
ing £, these orbits are deflected away from the triple 
collision along an unstable direction U^ Ze of the TCP 
fixed point. The fraction of orbits reaching the surface E 
among those leaving the chaotic scattering region scales 
thus as in (fTTj) now with exponent Xf^ e /|Ag°| . Similar 
mechanisms apply for the WR dynamics. 

The fraction A of two-electron trajectories making the 
transition from E back to E can thus in the limit E — > 0_ 
be estimated in terms of the stability exponents of the 
fixed points, also referred to as Siegel exponents 24] in 
celestial mechanics. One obtains A oc R^ 1 with fx as in 
(151) which can be written as 
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for the actual values of the stability exponents A, see 
7, 10]. The mean amplitude of the fluctuations in the 
quantum signal is thus asymptotically related to the frac- 
tion of phase space volume starting and ending at Rq. 

A numerical study of the full three-body quantum 
problem is still out of reach for energies E < In with 
iV ~ 15 [3,Q; we therefore chose a model system, namely 
eZe collinear helium first studied quantum mechanically 
in [25j. We calculated the cross section (TT|) directly in a 
large set of basis functions using the method of complex 
rotation and obtained a converged signal for N ~ 55. 
Semiclassically, we consider now the dynamics in the eZe 
space alone, which contains all the important parts re- 
garding the algebraic decay in the fluctuations. The 
number of basis functions used are scaled with energy 
to cover a fixed scaled region in R containing CTCOs 
with S/2ir < 20. Adopting the basis functions used by 
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FIG. 1: The Fourier spectrum of the fluctuating part of the 
eZe cross section rescaled according to (|12|l ; the circles denote 

— 1/2 

the position Sj and (relative) size of \M 12 \ 3 ' for CTCO's 
with S /2n < 13. Corresponding trajectories in configuration 
space are shown for the first 4 peaks. Inset: aji for N < 52. 

Puttner et al [l[ leading to a strongly banded Hamilto- 
nian matrix, it is possible to increase the basis size to 10 6 . 
Starting with an odd initial state 4>i , we obtain the cross 
section for the even parity eZe spectrum; its fluctuating 
part after numerically subtracting a smooth background 
is shown in Fig. [JJ The numerical value of the exponent 
fi is determined by rescaling the signal according to 

F(z) = \E\-»<j f i(z)/fuj (12) 

and testing the stationarity of the Fourier transform of 



F(z) in different energy windows (23[. The best value 
thus obtained is fi = 1.306±0.035 in good agreement with 
the theoretical prediction (JU) . (Note that the real parts of 
the exponents for 3-dim. helium and for eZe helium coin- 
cide as K/itur = 0). Furthermore, the peaks in the Fourier 
transform can be associated one-by-one with CTCO's in 
the eZe system, see Fig. [TJ We do not observe peaks as- 
sociated with the concatenation of different CTCO's or 
repetitions of single CTCO's. This is consistent with the 
expected suppression of orbit contributions traversing £ 
more than once as discussed earlier. Furthermore, we cal- 
culated the geometrical contribution to the coefficients aj 
in ([9]) directly from the matrix-elements M\2 by scaling 
out the leading order divergence; a clear correlation with 
the peak heights can be seen in Fig. [TJ Quantum contri- 
butions from the core region are thus indeed roughly the 
same for all CTCOs. 

In conclusion, we show that the fluctuations in the to- 
tal photo-ionisation cross section below the double ion- 
isation threshold follow an algebraic law with a novel 
exponent which can be written in terms of stability ex- 
ponents of the triple collision. Our findings are verified 
numerically for a collinear model systems; we further- 
more predict that the algebraic decay law is valid for 
the physically relevant 3 dimensional cases with an addi- 
tional contribution from the WR dynamics for Z > 9/4. 
Our findings will provide new impetus for experimental- 
ists and theoreticians alike to study highly doubly excited 
states in two-electron atoms. 
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